Localization and entanglement through scattering measurements 
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We study the emergence of locahzation and entanglement in many-body systems as a result 
of scattering measurements. We show that consecutive scattering measurements on a many-body 
system can produce superposition states in position space. The resulting states are signaled in the 
scattering distribution and these states are robust under continued scattering. We also show that 
the average scattering distribution for an ensemble of experiments gives the initial state scattering 
distribution even when the initial state is completely destroyed by the scattering process. 

PACS numbers: 03.65.Ta 37.10.Jk 03.65.Nk 
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Scattering of particles from matter is a fundamental 
method used to reveal a system's properties. A promi- 
nent example is Rutherford's revolutionary discovery of 
the atomic nucleus by scattering alpha particles from 
gold foil Scattering of neutrons has been used to 
probe the nature of exotic forms of matter such as su- 
pcrfluid helium and high temperature superconductors 
[3] , while scattering of x-ray photons has led to numer- 
ous discoveries about crystalline structure, such as the 
structure of DNA and scattering of optical photons 
has for example been used to probe correlations within 
Bose-Einstein condensates Neutral probe particles 
scatter from matter either as a result of density fluctua- 
tions in the system or by creating excitation in the sys- 
tem. Measuring the scattering pattern then yields valu- 
able information about the system's spatial structure and 
excitations, and furthermore interference between differ- 
ent scattering channels reveals information about corre- 
lations within the system. 

Typically in experiments and theoretical discussions 
involving scattering the effect of the probe on the sys- 
tem is assumed to be weak and the scattering pattern is 
calculated in the Born approximation from the initial sys- 
tem state [Ij. However during each scattering event en- 
ergy and momentum may be transferred from the probe 
to the system. Furthermore the act of measuring the 
probe's scattering pattern may produce quantum back- 
action on the state of the system. These effects can com- 
bine to significantly alter the state of the system, indeed 
for two particles, scattering of light leads to localization 
of relative position while the particles' momentum states 
become entangled 0] ■ 

Here we model the evolution of many-body states 
caused by consecutive scattering and measurement 
events, and show that for a system of multiple parti- 
cles this evolution can lead to many-body superpositions 
in position space. Such states may be useful for making 
measurements with enhanced precision and for investi- 
gating the boundary between the classical and the quan- 
tum. Our model also shows that the average scattering 
pattern for an ensemble of experiments will revert to the 
pattern predicted for the initial state even if this initial 
state is completely destroyed in each individual scatter- 
ing experiment. Experimental access to the initial state 
of system is hence possible even for systems that are eas- 



ily perturbed. 

We begin by examining the scattering problem from 
a general perspective. We consider an atomic system of 
identical particles with associated field creation and an- 
nihilation operators ^'''(r) and ^(r). The probe particles 
are assumed to have a definite momentum with associ- 
ated wave- vector kg . Scattering of the probe by the sys- 
tem is then associated with a momentum transfer to the 
system of Kk = S(ko — k/), where kj is the probe wave- 
vector after scattering. Measurement of the final probe 
wave- vector, for example by imaging in the far field, then 
transforms the system state according to 
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For an initial A^-particle state I ^/i) = / (iR(/)(R)|R), where 
R = {ri,...,rjv} gives the coordinates of the N parti- 
cles, the associated probability of observing a scattered 
probe wave- vector of k/ 
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The factor g is determined by the strength of the probe- 
system interaction and is assumed to be small so that 
multiple scattering of an individual probe particle is neg- 
ligible. For simplicity we assume g is independent of 
kj, which is the case when energy conservation restricts 
scattering so that |kf| ~ |ko| and the probe-system in- 
teraction is isotropic [5|. Anisotropic scattering, such as 
the dipole pattern in light scattering, can be taken into 
account by straightforward generalization. 

The probability of scattering to any final wave- vector 
is highest when all the particles are localized to have 
zero relative position and the scattering probability be- 
comes uniform. More generally the probability of scat- 
tering increases as the particles' relative positions be- 
come increasingly well localized, and observing a scat- 
tering event projects the state towards greater localiza- 
tion of the particles. The scattering distribution depends 

'^j^i e*'' ''j , which is completely determined by the 

relative positions of pairs of particles. The information 
gained from a measurement of the scattering distribution 
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is therefore limited to relative position and the measure- 
ment will then preserve superpositions of states that have 
the same set of relative position vectors. Each state in 
such a superposition will give the same scattering pat- 
tern and consecutive scattering measurements can push 
the system towards this type of superposition. A similar 
effect is observed in arrays of Bosc-Einstcin condensates, 
where measurement of atoms coupled out of the conden- 
sates results in superposition of relative phase . There, 
the relative phase localization was intrinsically limited by 
the number of atoms in the condensates, whereas here the 
number of scattering events and the resulting localization 
is unlimited. 

Apart from scattering we must also consider the situ- 
ation when the probe does not scatter from the system. 
Rather than leaving the system in the same state, the 
measurement of a non-scattering event also leads to a 
change in the state @ . This is because non-scattering is 
more likely from certain particle distributions and detect- 
ing a non-scattering event projects the wavefunction of 
the system towards these configurations. For a position 
eigenstate |R) the probability of non-scattering is deter- 
mined by the conservation of probe particles, that is the 
probe particles must either be scattered or not scattered, 
giving 



Pws(R) = |A(R)I' = 1 - 1^ y" dn^, 
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where the integration is over all scattering angles. De- 
tecting a non-scattering event then transforms the system 
state according to 
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and occurs with total probability / dR| A(R)0(R)|2. 
This projection favors basis states that have a lower prob- 
ability of scattering and hence works in the opposite di- 
rection to the projection following a scattering event. As 
we will sec this prevents the system from always ending 
in states with all the atoms localized with zero relative 
position. 

Consecutive scattering and non-scattering events lead 
to a dynamic evolution of the many-body state. The dy- 
namic scattering process can be simulated by the follow- 
ing procedure. Taking the initial state, we calculate the 
probability of non-scattering. A pseudo-random num- 
ber is then generated to determine if non-scattering oc- 
curs. If it docs then the projection in Eq. Q is applied. 
If instead scattering occurs then the random number is 
used to determine the angle of scattering according to 
the probability density in Eq. ((2), followed by projection 
of the state using Eq. ([T]). In either case the many-body 
state is then normalized and becomes the input state and 
the process repeats. 

We now examine an example simulation that displays 
the key characteristics of the dynamic scattering pro- 



cess. We consider matter- wave scattering from a one- 
dimensional lattice of bosonic atoms in a sinusoidal po- 
tential 0, which can be achieved in an optical lattice 
[l0| . We take the lattice to have M sites and to be ori- 
ented along the y-axis, while the initial wave-vector of 
the probe is in the x-direction. For simplicity we only 
consider scattering within the xy-plane, where the full 
three dimensional scattering is a straightforward gener- 
alization. 

At low temperature the atoms all reside in the low- 
est band of the lattice and the energy of the probing 
matter-waves can be arranged so that excitation of the 
atoms in the lattice to higher bands is negligible [1]. 
The atomic field operator can then be reduced to an 
expansion in terms of lowest band Wannier functions 
^{^) = '^jU!{r — Tj), where rj is the position of the 
jth lattice site We assume that the lattice poten- 

tial is strong enough that the overlap between neigh- 
boring Wannier functions is negligible. The state of 
the system can be expressed in terms of a number ba- 
sis = \{n^^\j = where nj"^ is the 
number of atoms at site j and u uniquely identifies each 
basis state. After the mth scattering event we expand 
the state as l^*™) = J2u Cu"' !{"}«). Scattering from this 
state then occurs at angle 9 to the x-axis with probability 
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where I{e) ^ J dreM'^y''\w{r)\^ and k{9) = A:o(l - 
cos 6*, — sinS). 

Following a detection at 6 the many-body state is pro- 
jected into the new state 

|vI/„+i) = -^^cMJ^e^-'-WnJ-^IMO, (6) 



where Af normalizes the state. We note that the number 
basis states are eigenstates of this projection, and the 
scattering process will preserve any state that begins in 
a basis state. Moreover, some of the basis states produce 
the same scattering pattern as the relative positions of 
pairs of atoms in the lattice is the same, for example 
in the TV = M = 3 case, |201) and |102) result in the 
same light scattering. Superposition of these states are 
partially preserved by the projection, in that the weight 
of each state in the superposition remains the same after 
scattering but the phase relationship is changed. For 
initial states that are the superposition of all the basis 
states, the projection gives higher weight to the basis 
states where the relative position is zero for each pair 
of particles, that is the state where Ug = TV, for some 
lattice site s, and all other sites have zero occupancy. 
These states have the highest probability of scattering at 
any angle and detecting a scattering event makes it more 
probable that the atomic system is in one of these states. 
A sequence of scattering events with 9 then leads to 
a superposition of these states in the limit to — > oo. 
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FIG. 1. Development of an atomic state caused by 3000 scat- 
tering events for a lattice with three sites and three atoms, 
(a) Modulus squared of the basis coefficients of the state 
l^'m) • (b) Detected events, (c) Overlap of the many-body 
state with the initial ground state |^o). Parameters used 
are U/J = 0.05, gN = 0.1 and ko = it /a, where a is lat- 
tice site separation and U and J are the parameters of the 
Bose-Hubbard model. 



In most cases, however, non-scattering events prevent 
the system from ending with all particles on one site, 
where these occur with probability 
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Detecting a non-scattering event projects the state into 
the new state 



*™+i) = -^^c('")A„|{n}„), 
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leading the state towards configurations that produce less 
uniform scattering distributions. 

The evolution of the lattice many-body state can now 
be simulated using the above probabilities and projec- 
tions. To examine the dynamic process we look at the 
simple case of a three site lattice containing three atoms. 
In this case there are only ten basis states making it 
straightforward to track the development of the many- 
body state. We solve the Bose-Hubbard model numeri- 
cally to find the initial many-body state of atoms in the 
lattice [m . In Figure [1] we show one realization of the 
dynamic scattering process for the three site lattice in the 
superfluid regime of the Bose-Hubbard model. For this 
example we have set the coupling constant gN = 0.1 and 
as a result the vast majority of detection events result 
from non-scattering, however it only takes a few scatter- 
ing events for the system to start favoring a particular 
final state. In this example the many-body state pro- 
gresses toward a superposition of the states |201) and 
1 102), two states which produce the same scattering pat- 
tern. As discussed above, continued scattering from this 
end superposition does not change the constituent basis 
vectors but does change the phase of the superposition. 
In Figure [IJc) we see that detection at non-zero angle 
quickly reduces the overlap of the many-body state with 
the original ground state. We see that the overlap makes 
quantum jumps when a scattering event occurs and grad- 
ually evolves due to non-scattering events. 

By repeated simulation of the dynamic scattering pro- 
cess, we find that the simulations all settle into an end 
state after a number of scattering events. The end states 
are always superpositions of eigenstates of the scatter- 
ing projection that produce the same scattering pattern. 
For the three site, three atom case there are four types of 
end state: (1) superpositions of |300), |030), and |003), 
(2) superpositions of |210), |120), |012) and |021), (3) 
superpositions of |201), and |102), and (4) |111). Each 
superposition is associated with a different scattering dis- 
tribution, which could be used to distinguish the super- 
positions in an experiment. 

In the Mott insulator regime of the Bose-Hubbard 
model the initial ground state is dominated by the |111) 
basis state and the most likely outcome of the dynamic 
process is to drive the state completely into the |111) 
state. In fact we find that the probability the end state 
of the dynamical scattering process is |111) is equal to 
|(lll|^'o)P, where |\l/o) is the initial state. This general- 
izes for final states that are superpositions, for example 
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FIG. 2. (a) Proportion of 10000 simulations that end in a 
final state superposition of (1) |300), |030), and |003), (2) 
|210>, |120>, |012> and |021>, (3) |201), and |102), and (4) 
I 111) at various values of U/J. The initial probabilities of 
finding these basis states in the initial ground state as shown 
by the diamonds, (b) (points) Number of event detections 
Nd{0) after 10000 simulations of 1000 detection events each, 
where the angular range is divided into 600 bins for event 
counting, (lines) Scattering distributions predicted from the 
initial ground states. Parameters used are gN = 0.5 and 
ko = ir/a. 



the end state superposition of basis states |201) and |102) 
occurs with probabihty | (201|«'o) P + |(102|4'o)p. In Fig- 
ure [H^a) we plot the proportion with which each possible 
final state superposition is represented in an ensemble of 
10000 simulations for various ratios of the on-site inter- 
action strength U and tunneling energy J. We see that 
the proportions match the initial probabilities of finding 
the relevant basis states in the initial ground state. The 
Mott state almost always ends in the |111) state, while 



for the superfluid, the end states are distributed across 
all possible outcomes. These results are independent of 
the coupling parameter g, which only affects the rate at 
which the state approaches an end state. 

Because the final state proportions are the same as 
the initial basis state probabilities, we expect that the 
scattering intensity distribution will equal that given by 
the initial ground state scattering probability distribu- 
tion if we average over an ensemble of scattering exper- 
iments. This occurs even though the initial state can be 
completely changed in the scattering process. In Figure 
[SJb) we show the intensity distributions for two systems, 
one in the superfluid regime and the other in the Mott 
regime, for an ensemble of 10000 scattering simulations 
with 1000 detection events each. The scattering inten- 
sity distributions for each ensemble match the scattering 
probability for the initial ground states, despite that in 
the superfluid case the individual realizations do not give 
the correct distribution. For light scattering from atoms 
in a lattice we expect a similar result, however in this case 
contributions from higher bands must be accounted for 
[l3| . Light scattering from optical lattices into a cavity 
has also been suggested as a method for creating partic- 
ular many-body states [l^ |. 

In this treatment we have not included the effect of 
time development due to the atomic Hamiltonian, un- 
der the assumption that the scattering process occurs 
on a time scale short compared to the tunneling and in- 
teraction time scales. For slower scattering this could 
be taken into account by evolving the system using the 
atomic Hamiltonian in between the scattering events. 

Our example of scattering from atoms in a lat- 
tice shows how scattering and measurement pushes the 
atomic many-body state towards eigenstates of the scat- 
tering projection. In some cases this leads to super- 
positions of many-body states in position space, where 
the constituent states give the same scattering pattern. 
These states are robust under continued scattering, un- 
like other more fragile superpositions where interactions 
with the environment lead to decoherence [l5| . The con- 
stituents of the superposition are revealed by the scatter- 
ing pattern that results after multiple scattering events, 
and scattering could be used to generate these entangled 
states in a non-deterministic fashion. It is also apparent 
that scattering distributions reveal the initial state scat- 
tering distribution, even when the initial state is signifi- 
cantly perturbed, provided the distributions are averaged 
over a number of experimental runs. 
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